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Abstract 

A new Lagrangian formulation with complex currents is developed and 
yields a direct and simple method for modeling three-phase permanent- 
magnet and induction machines. The Lagrangian is the sum a mechan- 
ical one and of a magnetic one. This magnetic Lagrangian is expressed 
in terms of rotor angle, complex stator and rotor currents. A complexifi- 
cation procedure widely used in quantum electrodynamic is applied here 
in order to derive the Euler-Lagrange equations with complex stator and 
rotor currents. Such complexification process avoids the usual separation 
into real and imaginary parts and simplifies notably the calculations. Via 
simple modifications of such magnetic Lagrangians we derive new dynam- 
ical models describing permanent-magnet machines with both saturation 
and saliency, and induction machines with both magnetic saturation and 
space harmonics. For each model we also provide its Hamiltonian thus its 
magnetic energy. This energy is also expressed with complex currents and 
can be directly used in Lyapunov and/or passivity based control. Further, 
we briefly investigate the observability of this class of Euler-Lagrange mod- 
els, in the so-called sensorless case when the measured output is the stator 
current and the load torque is constant but unknown. For all the dynam- 
ical models obtained via such variational principles, we prove that their 
linear tangent systems are unobservable around a one-dimensional family 
of steady-states attached to the same constant stator voltage and current. 
This negative result explains why sensorless control of three-phase electri- 
cal machines around zero stator frequency remains yet a difficult control 
problem. 
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1 Introduction 



Modeling electrical machines with magnetic-saturation and space-harmonics ef- 
fects is not a straightforward task and could lead to complicated developments 
when a detailed physical description is included (see, e.g., [31 [I]). Even if such 
effects are not dominant they could play an important role for sensorless con- 
trol (no rotor position or velocity sensor). For a permanent-magnet machine, 
the rotor position will be unobservable without saliency. For standard models 
of induction machines (no magnetic saturation, no space-harmonics) the rotor 
velocity is always unobservable at zero stator frequency [U H|. A global ob- 
servability analysis based on such standard models is given in [7J. In this note 
we develop a systematic method to include into such standard control models 
magnetic-saturation, saliency and/or space-harmonics effects. Our initial moti- 
vation was to see whether such non-observability is destroyed by such modelling 
changes or not. It appears that any physical-consistent model admits the same 
kind of observability deficiency at zero stator frequency. 

By physically consistent models we mean Lagrangian-based models. We pro- 
pose here an extension of Lagrangian modeling of three-phase machines with 
real variables (see, e.g. [9]) to complex electrical variables. It is directly inspired 
from quantum electro-dynamics where Lagrangian with complex generalized po- 
sitions and velocities are widely used (see, e.g., [4], page 87). We obtain, from 
such Lagrangian functions, physically consistent and synthetic Euler-Lagrange 
models directly expressed with complex stator and rotor currents. Such mod- 
eling method by-passes the usual detailed physical descriptions that are not 
easily accessible to the control community. Here we propose a much more direct 
way: it just consists in modifying the magnetic part of the Lagrangian directly 
expressed with complex currents and then in deriving the dynamic equations 
from the Euler-Lagrange with complex variables. We obtain automatically the 
dynamics of the electrical part as a set of complex differential equations. We 
suggest here simple Lagrangians modeling simultaneously magnetic-saturation, 
saliency and space-harmonic effects. The obtained dynamics extend directly the 
ones used in almost all control-theoretic papers and include also more elaborate 
ones that can be found in specialized books such as [1]. 

For permanent-magnet three-phase machines, the general structure of any 
physically consistent model including magnetic-saturation, saliency and other 
conservative effects is given by equations ((4]) with magnetic Lagrangian C m 
and Hamiltonian Ti m related by ([5]). For induction three-phase machines the 
physically consistent models are given by (|I5p where the magnetic Lagrangian 
C m is related to the magnetic energy Ti. m by (|f 6p . Such synthetic formulation of 
the dynamical equations is new and constitutes the first contribution of this note. 
We propose here natural modifications of the standard Lagrangian to include 
magnetic-saturation, saliency and space-harmonics, derive the corresponding 
dynamical equations and magnetic energies that could be used in the future to 
construct controlled Lyapunov functions and/or storage functions for passivity- 
based feedback laws. 

From a control theoretic point of view we just prove here that the severe ob- 
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servability difficulties encountered in sensorless control and well explained in [7] 
for the standard model resulting from the quadratic Lagrangian (| 13[) . remain 
present for models (|f 5[) where the magnetic Lagrangian is any function of the 
rotor angle, stator and rotor currents. Consequently, addition of magnetic sat- 
urations, saliency and harmonics effects, do not remove observability issues at 
zero stator frequency in the sensorless case (see proposition [I}. This observabil- 
ity obstruction has nether been stated for models with magnetic-saturation and 
space-harmonics of three-phase machines and constitutes the second contribu- 
tion of this notes. Contrarily to observability, non-observability is not a generic 
property and could be destroyed by generic and small changes in the equations. 
Since proposition [1] is based on the class of models derived from (fT]) or (TT5)) with 
arbitrary magnetic Lagrangian £ m , we prove here that any physically consis- 
tent model of three-phase machines where the non-conservative effects result 
only from voltage supply and Ohmic losses, such non-observability holds true 
around zero stator frequency. This means that non-observability around zero 
stator frequency is robust to generic and physically consistent modifications of 
the equations. As far as we know this negative and physically robust result 
is new. It indicates that sensorless control of three-phase electrical machines 
around zero stator frequency cannot be just addressed via refined physical mod- 
els but also requires advanced and nonlinear control techniques. 

In section[2]we recall the simplest model of a permanent- magnet machine and 
its Euler-Lagrange formulation based on the two scalar components of the com- 
plex stator current. Using the complexification procedure detailed in appendix, 
we show how to use complex representation of stator-current in Lagrangian 
formulation of the dynamics. This leads us to the general form of physically 
consistent models (equation |4|)). Finally we obtain, just by simple modifica- 
tions of the magnetic Lagrangian, physically consistent models with magnetic 
saturation and saliency effects (equation (flOf ). Section [3] deals with induction 
machines and admits the same progression as the previous one: we start with the 
usual (a, (3) model, describe its complex Lagrangian formulation, derive physi- 
cally consistent models (equations (|15[) ) and specialize them to saturation and 
space- harmonics effects (equation (f!H]) ). In section 21 we prove proposition [1] 
that states the main observability issues of these Euler-Lagrangian models at 
zero stator frequency. In conclusion we show how to transpose this modelling 
based on complex currents associated to a Lagrangian formulation to complex 
fluxes associated to a Hamiltionian formulation with complex generalized po- 
sitions and momentums. The appendix details the complexification procedure. 
It explains how to derive the Euler-Lagrange equations when some generalized 
positions and velocities are treated as complex quantities. Throughout the pa- 
per, we define models in (a, (3) frame, using the standard transformation from 
three phases frame (see, e.g., [5]). 
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2 Permanent-magnet three-phase machines 

2.1 The usual model and its magnetic energy 

In the (a,/?) frame (total power invariant transformation), the dynamic equa- 
tions read (see, e.g., [5]): 

/ _ ;\ „//- — \ 

(1) 




where 

• * stands for complex-conjugation, 3 means imaginary part, j = y— 1 and 
n p is the number of pairs of poles. 

• 8 is the rotor mechanical angle, J and tl are the inertia and load torque, 
respectively. 

• i s £ C is the stator current, u s £ C the stator voltage. 

• A = (Ld + L q )/2 with inductances Ld = L q > (no saliency here). 

• The stator flux is <f> a = Ai s + (f>e jnpe with the constant <j> > representing 
to the rotor flux due to permanent magnets. 

The Lagrangian associated to this system is the sum of the mechanical one C c 
and magnetic one C m defined as follows: 

£ c = ^ 2 , C m = ±\i a +te> n > e \ 2 (2) 

where i = <p/X > is the permanent magnetizing current. 

It is well known that (fT]) derives from a variational principle (see, e.g., [5]) 
and thus can be written as Euler-Lagrange equations with source terms cor- 
responding to energy exchange with the environment. Consider the additional 
complex variable q s £ C defined by = i s . Take the Lagrangian £ = C c + C m 
as a real function of the generalized coordinates q = (9, q sa , q s p) and generalized 

velocities q = {6,i sa ,i 8 p): 

C(q,q) = ^9 2 + - ((q aa +icosn p 9) 2 + (q sl3 + isinn p 6) 2 ) (3) 
Then the mechanical equation in (fT]) reads 

d fdC\ dC 



dt\d6 J 96 TL 

where — tl corresponds to the energy exchange through the mechanical load 
torque. Similarly, the real part of complex and electrical equation in (fT]) reads 

d ( dC\ dC 

dt{d q V ] -^ = Usa - Rsls 
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and its imaginary part 



d ( dC\ dC 

dt \dq s (3 ) dq s p 

since -j^- — -J^- — and q s — i s . The energy exchanges here are due to the 
power supply through the voltage u s and also to dissipation and irreversible 
phenomena due to stator resistance represented by the Ohm law —R s i s . 

2.2 Euler-Lagrange equation with complex current 

The drawback of such Lagrangian formulation is that we have to split into real 
and imaginary parts the generalized complex coordinates with q s = q sa + jq s p, 
(q sa and q s p real) and velocities q s — i s = q sa + jq s /3, (q sa and q s p real). We 
do not preserve the elegant formulation of the electrical part through complex 
variables and equations. 

Let us apply the complexification procedure detailed in appendix to the 
La,gra,ngi&n £(9 7 q sa ,q s p, 9, q sa ,q s p) defined in ([3|). The complexification process 
only focuses on q s and q s — i s by considering £ as a function of (9, q s , ql", 9, i s , i* s )'- 

Following the notations in appendix, n c — 1 with q c = q s , n r = 1 with q T = 9, 
S r = —tl and S c = u s — R s i s - Then according to (|2"5)) the usual equations |T]) 
read 

d fdC\ dC n d (dC\ 

dt{-^) = d9- TL > 2 dt{d?J =Us - Rsls 



sinC e |J = and J§ = f|. 

More generally, for any magnetic Lag reLligielli C^rn t licit is a real value function 
of 9, i s and 1* and that is ^ periodic versus 9, we get the general model (with 
saliency, saturation, space- harmonics, ...) of three-phase permanent-magnet 
machine: 

d / tA \ dC m „ d ( dC 

dt 



We recover the usual equation with (f> s = corresponding to the stator flux 

and to the conjugate momentum p c of q c as shown in appendix. According 
to (j2"5|) . the Hamiltonian H is the sum of two energy: H = H c + H m - The 
mechanical kinetic energy H c = ^9 2 and the magnetic energy 

H m {9,t s ,l*) = -q^—Is + ~Q^ % *s ~ An- (5) 

The standard model (fT]) derives from a magnetic Lagrangian of the form C m — 
■| \i s + ie jn p e \ with A and 1 are two positive parameters. Its corresponding 
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magnetic energy reads TL m = (|z s | 2 — z 2 ) . We recover the usual magnetic 
energy -||«s| 2 up to the constant magnetizing energy -|z 2 . 

In (|3|) , many other formulations of C m are possible and depend on partic- 
ular modeling issues. Usually, the dominant part of C m will be of the form 
^ \i s + ie 3ripe \ (A, i positive constants) to which correction terms that are 
"small" scalar functions of (9,i s ,i* s ) are added. 

2.3 Saliency models 

Adding to C m the correction — fSft (i 2 e _2:, ™ p9 ) with |^| < A means real part) 
provides a simple way to represent saliency phenomena while the dominant 
part of the magnetic Lagrangian (and thus of the dynamics) remains attached 
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to 4 Us + ie jnp \ . With a magnetic Lagrangian of the form 



- ^ I [i j ' 



''» e y + {i s e-^ 9 ) 2 ) (6) 



where A = [Ld + L q )/2 and /i = (L q — Ld)/2 (inductances Ld > and L q > 0), 
equations ^ become (At — <fi) 

^ (JO) = n p 3 ((A*: + ^e-^ 9 - ^ s e- 2 ^ e )i s ) - r L 

j ' (7) 

- (\i s + ^ lp9 - ^e 2 ™ 8 ) = u s - R s i s 
at 

and we recover the usual model with saliency effect. In this case the magnetic 
energy is given by: 

n m = ±(M 2 -?) - J((,: e ^) 2 + (, s e-^) 2 ). (8) 

2.4 Saturation and saliency models 

We can also take into account magnetic saturation effects, i.e., the fact that 
inductances depend on the currents. Let us assume hrst that only the induc- 
tances A and fi in ([6]) depend on the modulus p = \i s + ie? np6 \. The magnetic 
Lagrangian now reads 



Ml 



/MK+^I) ((, ; ^) a +(,.e-^) a ). (9) 



The dynamics is given by (g]) with such C m . Denote A' = With ^| 



>-r— — ,„ wr A' and 54 = i s+le \^ m A', we get the following dynamical model 
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with both saliency and magnetic-saturation effects: 



(10) 



with A = A+ J 2 ^ an< i = 2 A* ■ ^ is interesting to compute 

the magnetic energy ri m from general formula (J5]): 

H ro = - \t a \ - -i 

\i s +ie^ B \\' / fl , ?n 

+ J : ! — i w jn P e +z^e jn 
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A ! + A ! 2|.„+»^«| 



^ l ^"" 8 ) 2 + (j S C- Jnpe ) 2 ) • (11) 



Such magnetic energy formulae are not straightforward. They are not ob- 
tained by replacing A and [i in the standard magnetic energy 1(5)) by A and M 
respectively. 



3 Induction three-phase machines 

We will now proceed as for permanent-magnet machines. Let us recall first the 
usual dynamical equations of an induction machine with complex stator and 
rotor currents. They admit the following form 

' j t (j6) = n p 3 (L m t* r e^ n " e i s ) - t l 

< j (L m (tr + i s e-^ e ) +L fr l r ) = -R r l r (12) 

^ ^ ( Lm (* s + l r eJnp0 ) +Lf s i s ) =u s - R s i s 

where 

• n p is the number of pairs of poles, 9 is the rotor mechanical angle, J and 
tl are the inertia and load torque, respectively. 

t i r £ C is the rotor current (in the rotor frame, different from the (d, q) 
frame) , i s € C the stator current (in the stator frame, identical to the 
(a,/?) frame) and u s e C the stator voltage (in the stator frame). The 
stator and rotor resistances are R s > and R r > 0. 

• The inductances L m , L f r and L f s are positive parameters with Lf r , Lf s <C 
L m . 

• The stator (resp. rotor) flux is <j) s = L m (i s + i r e plpe ) + Lf s i s (resp. 
<\>r = L m (« r + i s e~ jn ? e ) + L fr i r ). 
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3.1 Euler-Lagrange equation with complex current 

With notations of appendix, n c = 2 with q° — (i r ,i s ), n r — 1 with q r = 9, 
S c = (—R r i r ,u s — R s i s ) and S r = —r L . The Lagrangian associated to (fT2"]) , 
expressed with complex currents i r and i s , reads: 



C(9,e,i r ,i* r ,i s ,i*) = J -9 2 



The first term ^9 represents the mechanical Lagrangian and the remaining 
sum the magnetic Lagrangian C m . The dynamics (fT2")) read 

d fdC\ dC 

n d fdC\ n n d fdC\ 



dt\d6 09 



The magnetic Lagrangian C m has the following form that coincides here with 
the magnetic energy H. m : 

C m = I f\i s + ve^ e | 2 + ^fW\ 2 + ^M 2 . (13) 

More generally physically consistent model should be obtained with a Lagrangian 
of the form 

Am = ~e 2 +£ m (9,i r ,i* r ,i s ,i* s ) (14) 

where C m is the magnetic Lagrangian expressed with the rotor angle and cur- 
rents. It is ^ periodic versus 9. Any physically admissible model of a three- 
phase induction machine reads 

4 t (W) (15) 

where the rotor and stator flux are given by 



di* ' di* s ' 

In general the magnetic energy does not coincide with C m . It is given by 
that yields: 

1-1 (ft * *\ — ^^~" m I OC-m * 

'T.jtj yu, V j V' *«> — Ir ~r V 

+ ^7 ls + ^- Cm - (16) 
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3.2 Saturation models 



A simple way to include saturation effects is to consider that the main induc- 
tances L m appearing in (T3)) depends on the modulus p = \i a + i r e> npe \. Thus 
we consider the following magnetic-saturation Lagrangian: 



r — 



L m (\i s +i r e^" e \) 



Since (L' m 



dL m \ dC m 
dp I dd 



dL r< 



-jn p 6 



+ 



LL 



dL,, 



^V< + ^V:. (17) 



di* r 2\i s + i r e^ e \ m ' di* s 2\i s +i r e^p e \ m ' 
the saturation model (formula (fT5|) with C m given by (fl7|) ') reads 
d 



dt 
d 



(A m (i r + i s e~^ e ) +L fr i r ) 



(18) 



^ — (A m (t s + ire^ 6 ) +L fs i s ) =u s - R s i a 



with A*, 



L m + 



dt 

i s +i r e 3 "P 



Jnp e\ 



-L' function of \i s + i r e plp6 1 . We recover here 



usual saturation models (see, e.g., pQ, page 428). Notice the similarity with 
permanent-magnet machines and (fTOf . Following (fl6|) , the associated magnetic 
energy reads then 



L, 



\i s + i r e^p e \ L' 



jn p U 



L 



Notice also that such similarity is not complete and could be misleading: con- 
trarily to (|18[) that can be derived intuitively form (|12p by replacing L m by A m , 
the above magnetic energy is not provided by (fl~3|) with L m replaced by A m . 



3.3 Space-harmonics with saturation 

We can take into account space-harmonic effects by adding their contribution 
to the magnetic Lagrangian (fT7|) . According to pQ, page 298, the iron path in 
general gets shorter as the harmonic order gets higher. Thus saturation effect 
has relatively smaller influence on the spatial harmonics. Following [5], the 
Lagrangian C v of harmonic v is 

Cv = y {i s i* r e-i< T » l,n *- + i\i r e^ m P B ) (19) 
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with L„ a small parameter (\L U \ <C L m ) and with a v = ±1 depending on 
arithmetic conditions on v (see [5], equations (25) to (29)). The total magnetic 
Lagrangian now reads 

Cm = (.. + Ve^ e ) + <e~^) 

+ + ^M s + ^ (i.*^-**""* 9 + z^ r ^^ 9 ) . (20) 

Now the saturation model (fT8|) is changed as follows: 

^ = n p 9 ((A m e-^ 9 + LvWfT*""™**} l* T l 3 ) - Tl 

j t (A m (l r + t s e-^ + Lf r Xr + IkI*-*""™* 8 } = -Rrl r (21) 

^A m + tre 3 "" 9 ^ + L fa t a + Lvire^""™" 9 ^ = U a — R s ls 
Following (fH)]) , the associated magnetic energy reads then 

Several space-harmonics can be included in a similar way. Moreover saturation 
of space-harmonics can be also tackled just by choosing L u as a function of 
\i a + i r e~ 3n t' e \. As far as we know such explicit models including magnetic 
saturation and space harmonics have never been given. 

4 Observability issues at zero stator frequency 

The sensorless control case is characterized by a load torque T£ constant but 
unknown, control inputs u s and measured outputs i s . Models derived from ([4J 
for permanent-magnet machines (resp. from (j 1 5|) for inductions machines) can 
be always written in state-space form 

j t X = f(X,U), Y = h(X) (22) 

where X = (t l ,0, 0, »(t a ),9f(t,)) (resp. X = (t l ,6, 6, U(i r ), $(i r ), R(t,), 9f(»,))) 
with U = QSt(u s ), S(w 8 )), Y — (5ft(z s ), 3(« s )) and = 0. A stationary regime 
at zero stator frequency corresponds then to a steady state (X, U, Y) of (|2"2"j) 
satisfying f(X, U) — and Y — h(X). The tangent linear system around this 
steady state is then 

^-x = Ax + Bu, y = Cx (23) 
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where A = §l(X,U), B = §fj{X,U) and C = §^(X). If we assume that 
the linearized system (l23|) is observable, the Kalman criteria implies that the 

rank of the matrix ( ^ ~\ must be equal to dim(X). If it is the case, the 

mapping X i— > (f(X,U),h(X)) is maximum rank around X. This maximum 
rank condition just means that the set of algebraic equations characterizing 
the steady-state from the knowledge of U and Y, f(X, U) = and h(X) = Y 
admits around X the maximum rank dim(X). Such rank is not changed by any 
invertible manipulations of this set of equations characterizing the steady-state 
from the knowledge of the input and output values, U and Y. Putting the 
implicit Euler-Lagrange equations ((H) and (|15p into their explicit state-space 
forms (|22p involves such invertible manipulations. 

For permanent-magnet machines described by ([2} , this set of equations yields 
the following mapping 

• 8C 

(t l ,6, 9, %.) ^(0,9,-^-r L ,U- R s i s ,i 8 ) 

where C m depends on 9 and i s . Its rank should be maximum, i.e., equal to 
5. This is not the case since its rank is obviously equal to 4. For induction 
machines described by (|15j) . the mapping is 

(r L ,0,9,i r> i s ) h-> (0,0, ~ T L ,-R r i r ,U - R 8 i s ,i s ) 

where C m depends on 9, i s and i r . Its rank is equal to 6 whereas the maximum 
rank is 7. The above arguments yield following proposition: 

Proposition 1. Any dynamical model of 'permanent-magnet machines ((H) (resp. 
induction machines (115|) ) is unobservable around zero stator frequency regime 
when the measured output is the stator current i s and the load torque is constant 
but unknown. By unobservable we mean that: 

• to any constant input and output u s and i s satisfying u s = R s t s correspond 
a one dimensional family of steady states parameterized by the scalar vari- 
able £ with 

~ t l — QQ* = ^, i s = l s for the permanent-magnet ma- 

chines, 

~ t l = ^§f-(S,,0,0,i s ,i*), 9 = £, i r = 0, t s — % s for induction ma- 
chines. 

• the linear tangent systems around such steady-states are not observable; 

5 Conclusion 

The models proposed in this note, (fTO]) for permanent-magnet machine and (f2"Tj) 
for induction drives, are based on variational principles and Lagrangian formula- 
tion of the dynamics. Such formulations are particularly efficient to preserve the 
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physical insight while maintaining a synthetic view without describing all the 
technological and material details (see [10] for an excellent and tutorial overview 
of variational principles in physics) . Extensions to network of machines and gen- 
erators connected via long lines can also be developed with similar variational 
principles and Euler-Lagrange equations with complex currents and voltages. 

In this note we have put the emphasis on currents and thus Lagrangian 
modelling. Since flux variables are conjugated to current variables, Hamilto- 
nian modelling is also possible when fluxes are used instead of currents. For 
permanent-magnet machines, the Hamiltonian counterpart of Lagrangian mod- 
els (Q| reads: 

d ( A dH m d 8H m 

J9 = t l , —(p s =u s - 2R S 



dt\ ) 89 " dr" * * <9</>* 

where the magnetic energy H m is considered as a function of the rotor angle 
9, the stator flux 4> s and its complex conjugate <fi*. The stator current i s cor- 
responds then to . For induction machines, the Hamiltonian formulation 
associated to (1151) becomes 



d , 0f? dH m d dH r , 



d ( T0\ dUm 



dV' ' debt ' dt' 



The magnetic energy H m now depends on 9, the rotor flux d) r and its com- 
plex conjugate (f>*, the stator flux (j) s and its complex conjugate </>*, The rotor 
(resp. stator) current is then given by 2 ^'" (resp. 2 ^" )■ As for Lagrangian 
modelling, one can modify the magnetic energies of the standard models |T]) 
and (|12[) to include, for example, magnetic-saturation or space-harmonics ef- 
fects. This yields new formula expressing H m as function of angle and fluxes. 
The corresponding flux-based models are then given by the above equations. 

[Lagrangian and Hamiltonian with complex variables] It is explained in [H 
page 87] how to use complex coordinates for Lagrangian and Hamiltonian sys- 
tems. Here we propose a straightforward extension where the complexification 
procedure is only partial. Such extension cannot be found directly in text-book. 
In the context of electrical drives, such complexification applies only on electrical 
quantities whereas mechanical ones remain untouched. 

Assume we have a Lagrangian system with generalized positions q € K™, 
n > 3 and an analytic Lagrangian C(q, q). Let us decompose q into two set of 
components: 

• the first set q c = (qi , . . . , (?2n<= ) with < 2n c < n will be identified with n c 
complex numbers q% = qzk-i + JQ2k, k = 1, . . . , n c ; 

• the second set q r — (q2n c +i, ■ ■ ■ , Qn) gathers the n r = n — 2n c components 
that will remain untouched and real. 
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Thus we can identified q with (q c ,q r ) where q c £ C™ and q r £ M. n . Since the 
Lagrangian C is a real- value and analytic function, it can be seen as an analytic 
function of the complex variables q c , q c , q c *, q c * and of the real variables q r 
and q r (q c * corresponds to the complex conjugate of q c ). This function with be 
denoted by C(q c , q c *, q r ', q c , q c * ,q r ) an d is equal to C(q, q) where 



91 + « 9l - ?1 fee + <7n<= <7n<= - 



2 ' 2j 1 1 2 ' 2j 
Let us consider the Euler-Lagrange equations 
d ( d£\ dC 



, Ql , • • • , 



*tWJ = ^ +s «- <24) 

where the -terms correspond to non conservative energy exchanges with the 
environment. Similarly to q, wc decompose S = (Si, . . . , S n ) into S c £ C" and 
g jjn . g j g identified with (S c , S r ). We will reformulate these equations with 
C and its partial derivatives. For k = 2n c + k' with k' = 1, . . . ,n r , they remain 
unchanged j| (^§^-^ — -§^r~ + ^X'- The two scalar equations corresponding to 
k = 2k' — 1 and k = 2k' with k' — 1, . . . , n c yield to the following single complex 
equation ^ (2-J^hj = 2-J=^- + ££, since we have the identities 

„ <9£ <9£ dC n dC dC dC 
2tt- = tt Jtt— , 2— = + j- 



dql, dq 2 k'-i dq 2k - dqff dq 2k '-i dq 2 k> 

and similarly 

„ dL dC dC n dC DC dC 

2ttt- = ttt J-—- 2— = — + J- 



dqt dq 2 k'-i dq 2 k> dq?t dq 2k > 



'42k' 



This provides the following complex formulation of the real Euler-Lagrange 
equations (|24|) 



,,n c (25) 



tA^f)=^-F+ S 1 k = l,...,n\ (26) 




dt ydq r k J dq k 

In the usual complexification procedure ([H page 87]) the coefficient 2 appearing 
in the above equations does not appear. This is due to our special choice 
Qk = 92k— l + Jl2k instead of the usual choice q k = q2k - 1 J_ n ' 2k , This special 
choice preserves the correspondence, commonly used in electrical engineering, 
between complex and real electrical quantities 
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Let us assume that, for each q, the mapping q i— > is a smooth bijection. 
Then the Hamiltonian formulation of (j24|) reads 

d dH d m _ , 1 

^7* = a - ' 17^ fe = _ « hAfc ' k = 1j.--,« ( 2 <) 

with 7Y = Q£r ■ q — £ and p = Let us decompose p into p c S C" and 
p r G M n ' . Then p r = J^p and p c = 2^p§r. Simple computations yield another 
derivation of the Hamiltonian from C directly: 

~ dC dC dC ~ 

H = if^q c + ^q c * + ^q r - £ (28) 

dq c dq c * dq r V ; 

where H denotes the Hamiltonian H when is a considered as a function of 
(q c ,q c *,q r ,p c ,p c *,p r ). Then (27]) becomes 

d c „ dH d c dH c c 

a! r <97i d r „ r , r , . 
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